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We calculated the contribution of Coulomb energy to the spacing between the ground and first 
excited state of -^^^Th nucleus as a function of the deformation parameter 5. We show that despite 
the fact that the odd particle is a neutron, the change in Coulomb energy AUc between these 
two states can reach several hundreds keV. This means that the effect of the variation of the fine 
00 . structure constant a = e^ /he may be enhanced /SUc /E ~ lO'' times in the E =7.6 eV "nuclear 

^^ ' clock" transition between the ground and first excited states in the '^^^Th nucleus. 
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I. INTRODUCTION 



Unification theories applied to cosmology suggest a possibility of variation of the fundamental constants in the 

expanding Universe (see e.g. review ^]). There are hints of variation of a and mq^^/ Kqcd in quasar absorption 

spectra, Big Bang nucleosynthesis and Oklo natural nuclear reactor data (see fSj and references therein). Here Aqcd 

is the quantum chromodynamics (QCD) scale, and rUq and We are the quark and electron masses. However, the 

majority of publications report only limits on possible variations (see e.g. [H, @) S S Q ) • A very sensitive method 

O ■ to study the variation in a laboratory consists of the comparison of different optical and microwave atomic clocks 

^ ' (see e.g. measurements in [3, [3, M, [l3, [ill, U2, [l3, H, [HI). An enhancement of the relative effect of a variation 

^, can be obtained in a transition between almost degenerate levels in Dy atom [lg|. These levels move in opposite 

directions if a varies. An experiment is currently underway to place limits on a variation using this transition [l7|, 

^sj ■ but unfortunately one of the levels has quite a large linewidth and this limits the accuracy. An enhancement of 1-5 

^ ' orders exists in narrow microwave molecular transitions [l8|. Some atomic transitions with enhanced sensitivity are 

00 . listed in Ref. [19| . In Ref. [20| it was suggested that there may be a five orders of magintude enhancement of the 

J~: variation effects in the low-energy transition between the ground and the first excited states in the ^^^Th nucleus. 

i^ [ The exsitence of the enhancement was confirmed in [21|. This transition in ^^^Th was suggested as a possible nuclear 

clock in Ref. [22] ■ Indeed, the transition is very narrow. The width of the excited state is estimated to be about 

t^ I 10"'' Hz [23( (the experimental limits on the width are given in [23l); The latest measurement of the transition 



energy [2^ gives 7.6 ± 0.5 eV, compared to earlier values of 5.5 ± 1 eV [231 ^^-'^ 3. 5 ± 1 eV 2^. This makes ^^^Th a 

possible reference for an optical clock of very high accuracy, and opens a new possibility for a laboratory search for 

the variation of the fundamental constants. Large interest in experiments exploring this possibility was expressed, 

^J^ ■ for example, in E. Hudson and E. Peik talks at special workshop (14-18 July 2008, Perimeter Institute) devoted to 

K^ , the variation of the fundamental constants and privite communications of E. Peik, P. Bciersdorfer, Zheng-Tian Lu, 

;_j D. DeMille, D. Habs and J. Torgerson. This motivated us to perform an independent calculation of the sensitivity of 

Cu the transition frequency v to the variation of a based on the following relation suggested in Ref. [29| : 

hSv = AUc— (1) 

a 

where AUc is the contribution of the Coulomb interaction to the energy spacing between these two levels. 

II. SHIFT OF NUCLEAR COULOMB ENERGY 

The calculation of the difference in total Coulomb energy between ground and excited nuclear states can hardly 
be done within existing nuclear models with the required accuracy. Instead, the calculation of the correction to the 
Coulomb energy when one proton or one neutron is added into a single particle orbital \v) is a much simpler problem. 

The ground state of ^^^Th nucleus is [NuzAJ^] = [633 5/2"'"]; i.e. the deformed oscillator quantum numbers are 
N — 6, Hz = 3, the projection of the valence neutron orbital angular momentum on the nuclear symmetry axis 
(internal z-axis) is A = 3, the spin projection S = —1/2, and the total angular momentum and the total angular 
momentum projection are J = U = A + T, = 5/2. The excited state is [Nuz-^- J^] = [631 3/2"*"]; i.e. it has the same 



N = 6 and n^ = 3. The values A = 1, S = 1/2 and J = 17 = 3/2 are different. With relation to the variation of alpha, 
we shall discuss below the contribution of the Coulomb interaction to the energy spacing between these two levels. 

We start from an even ^^^Th nucleus and calculate the shift of nuclear Coulomb energy when an odd particle is 
added into single particle state \i>). In an even nucleus the Coulomb energy can be presented as 

where Pc{r) is the charge density. For uniformly charged sphere we obtain the well known expression 

_3ZV 

which gives for Z=90, and A=228, the estimate, Uc ~1 GeV. The above equation can also be used to estimate the 
change in the Coulomb energy when one neutron is added. An addition of the neutron changes A, A— >A+1, therefore 

that gives for AUc the value of the order of -1 MeV. This estimate gives some average macroscopic value for AUc 
because in reality it depends on the state {v) where we put the neutron. 

Adding an odd particle to the state \iy) we obtain a change in the charge density 

{u\p,{r)\u)^p,{r)+Sp'"'{r), (3) 

where pc{r) is the charge density of a neighbor even nucleus. For the change in the Coulomb energy Eq.Q we obtain 

AU'c^ j d'rU,{r)5p'"'{r), (4) 

where Uc{r) is the single particle Coulomb potential 

c^cW = /rfV^. (5) 

Although ■^■^^Th is a deformed nucleus, for the estimates we shall neglect a small quadrupole component in the charge 
density Pc{i')- Effects of deformation will be discussed below. In approximation of uniformly charged sphere 

Uc{r) = { (6) 

[^ r>R, 

where i? is a Coulomb radius of ^^^Th, R — 1.2^^/^fm. We use this very simple form in order to obtain an estimate 
for the change in the Coulomb energy when going from one nuclear state to another state, and not in the absolute 
value. 

The ground state of the nucleus ^^^Th is polarized due to the addition of the nucleon to the ^^^Th core. The added 
nucleon introduces an admixture of the monopole state (which is a radial excitation) [30|, |3l| and the nuclear state is 
now: 



|OV) = v/l-e2|O+;0,)+e|M;0,), (7) 

where |0+) and \M) are the ground state and the isovector monopole in ^^^Th and 0^ is the wave function of the 
added nucleon. With 

Ef) — Em 

where F^ is a two-body nuclear interaction. 

Now we evaluate the Coulomb energy of this new polarized state |0V). The correction is: 

AU^ = 2e(0+; '^^l^ E ^cir, - r,)\M; 0,}, (9) 



Vc{r) is the Coulomb potential. This correction can be written in the form: 

AU^^2ejJsp{r)U,{r)£r. (10) 

Sp{r) is the transition density between the ground state and the monopole: 

The Z/A factor takes into account the fact that only protons interact via the Coulomb. 
The strong interaction matrix element we evaluate using a simple interaction: 

Fn = Fo<5(ri - ra). (11) 

Then: 

{0+;cJ),\Fn\M;cJ),)=Fq j 5p{r)\4>,{r)\''d\. (12) 

Finally: 

^^^, _ 2Z/AF, ^ 5p{r)Uc{r)<Pr j 5p{r)\Ur)\'<Pr 

The result obviously depends on 0^. 

One form used for the monopole transitions is: 



dr r^ dr 



6p^Ci3p + rJj^CY-y:jJ^, (14) 



where C is a constant and p is the ground state density. This form for the transition density corresponds to a change 
obtained by uniformly expanding the nucleus, but also can be derived from a sum rule [30l.l3ll|. This transition density 
has a node. Because of the node in dp the result for the polarization correction will be sensitive to the spatial form 
of the wave function cj)^, and even may change sign for different single-particle states. 

To obtain a quantitative estimate let us switch to the particle-hole basis \ph) instead of the single monopole 
resonance \M). The correction to the charge density introduced in Eq.Q can be presented as follows: 

<^P'"'W = E'^/^A5;'V'Ur)VA(r), (15) 

A, A' 

where ipxi'"') is the single-particle wave function, and Sp'^'^, is the correction to the proton density matrix of a neighbohr 
even nucleus, due to the odd nucleon added into the state \p). 

{u\al]alW)^{0\al]al\0) + SpTy. (16) 

Here (0|a^la^|0) is the proton density matrix of the neighbohr even nucleus. 
Substituting Eq. (fT5)l into Eq.(|4]) we obtain 

AC/^-1](C/c)va^Paa', (17) 

A, A' 

where {Uc)x'x ^^ ^ matrix element of the single-particle Coulomb potential Eq.®. 
The equation for (Jp^^v <^^^ ^^ found for example in [325 . It states 

SpZ' = S^xS.y + ^^f^ E {^h\FN\X2X')SpZx,- (18) 

Ai A2 

Here (AAi|i^jv|A2A') is the matrix element of effective residual interaction between quasi-particles. The first term in 
r.h.s. of Eq. p^ exists only for an odd proton. If the odd particle is a neutron, this term is absent. The second term 
is the polarization correction which exists both for an odd proton and an odd neutron. 



In the first order in proton neutron residual interaction we obtain for the correction to the proton density matrix 

n\ - nx 



6p 



XX' 



(AHF^kA'), 



where the state \v) refers to a neutron. The correction to Coulomb energy becomes 

nx - nx' 

V^cJx'X 
X,X' 

Eg. ((20)1 can be presented in the foUowing way, 



AU^C = J2(Uc)y> 



f-X - EA' 



-(AHF^kA'). 



(19) 



(20) 



AU!^^ j d\5Uc{r)\Mr)?- (21) 

Here we introduced the correction to potential energy SUc {r) that describes the reaction of the proton core to the 
presence of an additional neutron 



5Uc{r) - ^oE (f^^)A'A ^;^^^l(r)VA'(r-) 



A. A' 



CA - fA' 



(22) 



It follows from Eqs. (|20l22p that the correction SUc{r) does not depend on the neutron state. However, its mean value 
given by Eq. (PT|) does depend on the particular state \i'). The radial dependence of 6Uc{r) is shown in Fig.l. This 
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FIG. 1: Radial dependence of the correction to the Coulomb potential energy 

radial dependence shown is typical for a double charged layer. Due to small compressibility of a nucleus the change 
of a charge density can take place on nuclear surface only, and in good approximation it can be presented as a double 
charged layer. The mean value of the correction SUcir) over nuclear volume is negative and equal to -440 keV. Due 
to oscillating behavior, the expectation value of SUc{r) will obviously be state dependent. 

The calculations for ^^^Th were performed using Nilsson wave functions [33| for the valence neutron in order to 
account for deformation. In ^^^Th the ground state has the quantum numbers [633 5/2+] while the excited state has 
[631 3/2+]. Both of them belong to N=6 shell where the possible angular momenta in the Nilsson orbitals are L=2, 4, 
6. The shift of Coulomb energy calculated due to add of valence neutron into these states is shown in Table 1. Note, 



TABLE I: Shifts of the Coulomb energy due to odd neutron for N=6 shell orbitals. 



AUc (MeV) 



3d. 



3/2 



-0.741 



3d. 



■5/2 



-0.926 



2g7/2 



-0.848 



2g9/2 



-1.381 



l»ll/2 



+0.158 



li 



13/2 



+1.305 



that the shift for 1^13/2 and liii/2 becomes positive. The radial wave function for L=6 has a sharp peak near nuclear 
surface where 5Uc{r) is positive. The state liii/2 has smaller binding energy and its peak, compared to lii3/2 state, 
is shifted to larger r where 5Uc{r) becomes negative. This leads to a smaller value of the shift for liii/2 state. 

The difference between the first excited state with quantum numbers [631 3/2] and the ground state [633 5/2] in 
^^^Th is extremely small (< 10 eV) and can not be reproduced in Nilsson model at reasonable deformation 5 i2 0.2. 
With the standard parameters of the Nilsson model [S^] these two levels are crossed at the deformation 5 < 0.1. 




FIG. 2: Difference in Coulomb energies between the first excited and the ground states in ^^^Th as a function of deformation. 

Therefore, we calculated the difference in Coulomb energies ^Uj — AU^! at different deformations, as a function 
of the parameter 77 — 6/k, where k — 0.05 is the spin-orbit constant. The result is shown in Fig. 2 for < 77 < 6. 
This interval corresponds to deformations < 6 < 0.3. Note that the difference becomes negative at the deformations 
0.15 < S < 0.25. However, the position of this interval depends on the parameters of the Nilsson model, while 
the fact that the difference changes sign at some deformation remains in all sets of the parameters. The value of 
AUc as a function of deformation changes from 1.5 MeV at zero deformation down to -0.5 MeV a,t S = 0.3. In 
our approach we can hardly give a better estimate, due to simplified treatment of deformation. However, very small 
value of the Coulomb energy shift seems unprobable. Note that near the Nilsson level crossing the difference of the 
Coulomb energies is about 1 MeV. The shift of the order from several tens keV to several hundreds keV provides an 
enhancement of sensitivity to changes of a, as it was mentioned in [2d |. Indeed, the frequency shift from Eq.([T]) may 
be presented in the following form 

S^ ^^E^^^2A2x 10^^ Hz (23) 

lOOi^ey a ^ ' 

This shift is four to five orders of magnitude larger than the shift which was studied in optical atomic clocks [10, U^, 
llJj Hj llal ) • ^^ additional enhancement may be due to the very small width of this transition (10~^ Hz [23|). It is 
also instructive to present the relative enhancement: 

Su ^ AUc lOOKeV Sa^ . AUg Sa 

V IQQKeV 7.6eV a ' ^ lOOKeV a ^ ' 

These estimates confirm conclusion of Ref. [2(?| that experiments with ^^^Th have potential to improve sensitivity of 
laboratory measurements of the variation of the fundamental constants up to 8 orders of magnitude. 
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